Journal of Soft Computing in Civil Engineering 6-1 (2022) 88-109 


Journal of Soft Computing in Civil Engineering 


y oN Contents lists available at SCCE 
ww 


~ouYyan Journal homepage: www.jsoftcivil.com S C C @ 
Press 


Evaluation of Dimension Stone According to Resistance to Freeze— 
Thaw Cycling to Use in Cold Regions 


Reza Mikaeil'®, Akbar Esmaeilzade”®, Sina Shaffiee Haghshenas*®, Mohammad Ataei’®, 
es ee s s s e * 

Sina Hajizadehigdir’, Amir J afarpour’, Tae-Hyung Kim’, Zong Woo Geem® © 

1. Associate Professor, Department of Mining Engineering, Faculty of Environment, Urmia University of Technology, Urmia, 

Tran 

2. Assistant Professor, Department of Mining Engineering, Faculty of Environment, Urmia University of Technology, Urmia, 

Tran 

3. Ph.D. Candidate, Department of Civil Engineering, University of Calabria, 87036 Rende, Italy 

4. Professor, Department of Mining, Petroleum & Geophysics Engineering, Shahrood University of Technology, Shahrood, Iran 

5. M.S., Department of Civil Engineering, Faculty of Environment, Urmia University of Technology, Urmia, Iran 

6. Ph.D. Candidate, Department of Mining and Metallurgical Engineering, Yazd University, Yazd, Iran 

7. Professor, Department of Civil Engineering, Korea Maritime and Ocean University, Pusan 49112, Korea 

8. Associate Professor, College of IT Convergence, Gachon University, Seongnam 13120, Korea 

Corresponding author: geem@ gachon.ac.kr 


5) _https://doi.org/10.22115/SCCE.2022.325638.1398 


ARTICLE INFO ABSTRACT 

Article history: Freezing is one of the most effective natural and environmental factors 
Received: 20 January 2022 on the physical and mechanical characteristics of dimension stones. 
Accepted: 10 February 2022 Since, freezing is a destructive agent, thus causes the undesirable stone 


conditions and reduces quality and its efficiency. This study, it was 


aimed to evaluate and rank the dimension stones according to their 
changes in physical and mechanical properties due to freezing 
conditions. For this purpose, 14 rock types of the most widely used 
dimension stones in cold regions were collected and transferred to the 
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PROME HEE: laboratory to determine their physical and mechanical characteristics. 
PSO algorithm; 
Chastering In laboratory tests, standard samples of stones were prepared, and for 


all of the samples Uniaxial Compressive Strength (UCS), Durability 
Index (DD), Density (D), and Water absorption percentage (Wa) were 
determined before and after different freezing-thawing cycles. Then 
utility degree of studied stones in frost condition was assessed using 
the preference ranking organization method for enrichment of 
evaluations (PROMETHEE) multi-criteria decision-making method. 
The results of the study showed that samples of A3 (Piranshahr 
Granat), A10 (Hamadan_ black granite), A8 (Azarshahr yellow 
travertine), and A4 (Mahabad gray granite) are in order from the 
highest degree of desirability in a condition of freezing—-thawing and 
for use in cold climates are especially suitable for use in outdoor and 
urban spaces. In addition, the results of the laboratory were evaluated 
by the PSO algorithm for clustering analysis and com-pared with the 
ranking result by PROMETHEE. The results obtained demonstrated 
the proposed approach could be an efficient tool in the evaluation of 
the freezing phenomenon on physical and mechanical properties of 
dimension stones. 
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1. Introduction 


One of the most important factors that influences the efficiency of industrial activities is climate- 
related factors [1]. Temperature, as one of the climatic factors, has the greatest impact on mining 
activities, so that mining operations may be stopped in the event of severe temperature changes 
[2]. On the other hand, temperature changes cause the resistance of rocks to change [3]. Iran, 
located in Middle East, is famous for a diverse climate in the whole world that rapid temperature 
changes, the differences between hottest and coldest temp. During the day, hot sun and long 
freezes are sometimes diverse and conflicting climate properties of Iran [4]. The northwest and 
western provinces of Iran have mild summer and cold winters [5]. The high rainfall and cold 
climate are characteristics of the winter in these regions, and since the major potential and vast 
resources of decorative stones are located in these regions, stone owns a special place in the 
culture and art of housing in these regions. According to difficult climate conditions, facade 
stones suffer long-term freezing and multiple freezing—thawing cycles. This causes erosion and 
vast destruction of a variety of dimension stones of structures and reduces life and beauty [6,7]. 
A wide study on the effectiveness of freezing and thawing at the international level on the 
physical and mechanical properties of stones has taken place [8,9]. 


Many research works have been accomplished for the study of the impacts of freeze-thaw cycles 
on the physical and mechanical characteristics of dimension stones. Matsuoka (1990) studied the 
effectiveness of the freezing—thawing cycle according to water content in the stone [10]. In this 
regard, Nicholson and Nicholson (2000) inspected on decline mode for ten types of sedimentary 
rocks according to recrudescent cycles of freeze-thaw process [11]. In other studies, Mutluturk et 
al. (2004) offered a model to study the decay of stone during freezing—thawing cycles [12]. Then, 
Altindag et al. (2004) inspected it for the entirety damage of ignimbrite and specified the impacts 
of freeze-thaw cycles on Isparta ignimbrite, and this model was evaluated by these 
researchers[13]. Chen et al. (2004) study the welded porous of tuff from the viewpoint of the 
effect of water saturation to freeze-thaw cycles. The results indicate that both porosity and rock 
degrading significantly enhance when the saturation degree exceeds 70% [14]. Karaca (2010) 
inspected experimentally the effect of freeze-thaw cycles on the abrasion values [15]. Tan et al. 
(2011) studied the effect of freezing—thawing cycles on stone mechanical characteristics [16]. 
Also, Bayram (2012) presented a novel statistical model for estimating the percentage damage 
value in the uniaxial compressive strength of nine limestone examples after a freeze-thaw test 
[17]. Vehbi Gék¢e et al. (2016) reviewed the impact of the number of freezing—thawing cycles on 
physical and mechanical travertine characteristics using statistical studies made in Turkey [18]. 
The research was done in Iran in regard of effect of freezing and thawing on physical and 
mechanical properties of dimension stones, studies of Jamshidi et al. (2013), and also researches 
of Khanlari et al. (2015) in Iran, can be noted [19,20]. Cyclic freeze-thaw event in the cold and 
wet areas, change mechanical and physical characteristics of rocks. When water freezes, its 
volume increases up to 9%, and consequently by erecting dividing force on the sidewall of 
fractures and pores, increases pore pressure [21]. New discontinuities create or old fissure 
developed when pressure passes the level of the tensile strength of the rock pore; as a result, 
these changes can affect the stone’s mechanical strength and thus its durability [22,23]. Also, 
Jamshidi et al. (2016) presented a new physicomechanical index for predicting the mechanical 
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strength of travertine after a freeze-thaw test. The results of this research indicated that the 
physicomechanical parameter is well precise for predicting the mechanical strengths of travertine 
such as, uniaxial compressive strength and point load strength and Brazilian tensile strength after 
a freeze-thaw test, and consequently present a rapid durability evaluation [24]. In the other 
studies, Li et al. (2018) examined the time-dependent behavior of quartz sandstone and quartzite 
according to chemical degradation and freeze-thaw cycles [8]. Liping et al. (2019) were studied 
the physical index change and triaxial compression test of intact hard rock subjected to freeze- 
thaw cycles [25]. In another study, Peng et al. (2020), were evaluated the effect of rock freeze- 
thaw in cold-region tunnels [26]. Also, a deformation property of coarse-grained sulfate saline 
soil based on the freeze-thaw-precipitation cycle was studied by Zhang et al. (2020) [27]. 


In all these studies, researchers aimed to evaluate and research the coherence between different 
cycles of freezing and thawing and physical and mechanical properties of dimension stones 
which has resulted in valuable results. This study aimed in addition to evaluate the effect of 
quantity of freezing on most widely used dimension stones of the country, degree of desirability 
of studied stones using the PROMETHEE method. So that the most desirable dimension stone 
for use in urban structures is determined and provides a clear vision about its impact on the stone 
selection. It must be noted that there has not been much research on the ranking of dimension 
stones in different freezing—thawing cycles. Therefore, the necessity of research in this regard is 
undeniable. It should be noted that extensive studies have been done on the use of the 
PROMETHEE method in evaluating various mining engineering problems, some of which are 
mentioned in Table 1. 


Table 1 
Some of the studied using the PROMETHEE method [28,29,38—44,30-37]. 
Researcher(s) Year Description 
Dagdeviren 2008 Equipment selection 
Athawale & Chakraborty 2010 Facility location selection 
Prvulovic et al. 2011 Choice of systems for drying paltry-seeds and powder materials 
Bogdanovic et al. 2012 Mining method selection 
Abedi et al. 2012 Copper exploration 
Shiriskar & Patil 2013 Optimization of energy charges 
Balali & Zahraie 2014 Selection of building structural system 
Mikaeil et al. 2015 Sawability of dimension stone 
Mladineo et al. 2016 Priority setting in management of mine action projects 
Balusa & Singam 2018 Underground mining method selection 
Mikaeil et al. 2018 Geotechnical risks assessment in tunneling projects 
Ghadernejad et al. 2019 Selection of optimal coal seam 
Iphar & Alpay 2019 Underground mining method selection 
Sitorus et al. 2019 Choice problem in mining and mineral processing 
Mikaeil et al. 2020 Evaluating the coal seam mechanization 
Dayo-Olupona et al. 2020 Selection of emerging technology in surface mines 


Dachowski & Gatek 2020 Selection of the Best Method for Underpinning Foundations 
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According to the results of research conducted with the PROMETHEE method (Table 1), it can 
be said that this method can evaluate and rank in solving various problems of mining 
engineering. On the other hand, the results of the research show that this method is compatible 
with real-world conditions and provides reliable results for researchers and designers of mines. 


2. Field and laboratory studies 


This research contained two main parts as field and laboratory studies. In the field study, rock 
samples are provided from studied quarries located in Iran. In the laboratory study, the 
characteristics of physical and mechanical are determined before and after freeze-thaw cycling. 
Finally, the desirability degree of studied dimension stones is deter-mined using the 
PROMETHEE method and clustering technique. The flowchart of the re-search procedure is 
illustrated in Figure 1. 


eee eee of rock samples for laboratory Providing rock samples from studied 
eee eee quarries 


Transferring the samples to the 
rock mechanic lab 


Fig. 1. Flowchart of the study. 
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In this study, the numbers of dimension stones of west and northwest regions of Iran are checked. 
Frequency, locality, beauty, color, and price for these stones are the reason. So, 14 dimension 
stones are selected. The studied dimension stones are given in Table 2. The location of the 
studied quarries is shown in Figure 2. 


Table 2 
Studied dimension stones. 
Sample ID Rock-type City Province 

A, Cream travertine Khalkhal Ardabil 
Ad Walnut travertine Azarshahr East-Azerbaijan 
A3 Olive granite Piranshahr West-Azerbaijan 
Ag Gray granite Mahabad West-Azerbaijan 
As Red granite Naein Isfahan 
Ag Spring White granite Zanjan Zanjan 
Ay White travertine Tekab West-Azerbaijan 
Ag Yellow travertine Azarshahr East-Azerbaijan 
Ao Red marble Azarshahr East-Azerbaiyan 
Ato Black granite Hamadan Hamadan 
An Chocolate granite Khoramdare Zanjan 
An Cream marble Harsin Kermanshah 
A Cream Pink travertine Anarak Chopanan Isfahan 
Awa Pink marble Anarak Haftoman Isfahan 


Fig. 2. Flowchart of the study. 


S 


According to the nature of the freezing experiment, the need to control the uniformity of samples 
is of utmost importance. So, in this research to ensure the accuracy of the samples and to obtain 
uniform samples, instead of gathering samples from quarries, samples from stones factories in 


R. Mikaeil et al./ Journal of Soft Computing in Civil Engineering 6-1 (2022) 88-109 93 


Tabriz city were collected and transported to the rock mechanic laboratory. In total, a block with 
an average size of 40x40x50 cm? was selected from each type of stone. To begin studies in the 
laboratory, the 8 NX diameter cores and a length to diameter ratio of 5/2 were prepared. In all 
experiments, all steps have been according to Laboratory ISRM Standards [45]. According to 
ISRM, freezing experiment is performed at a temperature of —15°C, however in this study, 
considering that at least 20 nights a year at some of the cities in West and northwest of the 
country temperature reach —20°C even less, therefore, the temperature suitable for the simulation 
and realistic study of stones, laboratories were asked to test freezing and thawing in 5 courses at 
a temperature of —20°C. All tests were done once on the control samples and the samples which 
were spent freezing test and a week normally lost their moisture. It should be noted that the five 
of experiments for each sample before and after the freezing—thawing were conducted and the 
average of the results for each sample are shown in Table 3. Also, the results of tests on stone 
samples before and after the freezing—thawing are given in and Figures 3 to 6. 


Table 3 
Physical and mechanical properties of the studied stones before and after the freezing-thawing. 
D (gr/em’) Wa (%) UCS (MPa) DI (%) 


Sample ID Before After Before After Before After Before After 


freezing freezing freezing freezing freezing freezing freezing freezing 


Al 2.46 2.38 2.1 2.8 50.5 43.5 97.69 96.12 
Ao 2:52 2.48 2.72 2.95 48.5 41 98.87 97.67 
A3 3.18 3.13 0.41 0.43 105.5 95.5 99.38 98.68 
Ag 3.23 2.98 0.30 0.31 189 182.5 99.57 98.07 
As 3.31 2.87 1.02 1.1 1775 169 99.72 98.86 
Ag 2.87 2.69 0.96 1.05 121 111.5 99.43 98.49 
Ay 3.06 3 1.56 1.81 89 82.5 98.73 97.64 
Ag 2.56 2.31 2.61 2.76 33 32 98.43 97.29 
Ag 2.6 2.43 2.54 2.63 52.5 48 99.12 98.18 
Ato 3.26 3.01 0.25 0.265 173 166.5 99.56 98.88 
Au 2.95 2.74 0.98 1.05 133 126.3 99.31 98.1 
An 2.86 2.6 1.57 1.7 71.5 62 98.79 97.94 
Aj 3.07 2.72 1.94 2.03 73 65.5 99.14 98.23 
Aw 3.23 3.094 1.2 1.305 74.5 68 99.42 97.93 


D: Density, Wa: Water absorption percentage, UCS: Uniaxial compressive strength, DI. 
Durability index 
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As can be seen in Table 3, freezing—thawing more seriously affects the mechanical and physical 
properties of stones. Two important physical properties such as density and water absorption 
properties seem to be more effective by freezing—thawing. When the stones are wet, water 
permeates in stone pores and micro-cracks and causes freezing and the volume increase of 
existing water in these spaces increases internal tensions and affects the macroscopic structure of 
the stone. With ice thawing in these spaces, the structure of the stone releases from ice pressure, 
and very mild tension relief is created in the stone. This phenomenon increases stone micro- 
cracks and gets flaky inside and outside. A rapid freezing—thawing cycle toughens up the 
destruction of freezing. All the process alters the physical properties of stone, especially the two 
mentioned parameters and decreases in mechanical properties including compressive strength 
and durability index. 


Before freezing After freezing 


3.4 
3.2 
3 
& 2.8 
ob 
~ 2.6 
a 
2.4 
2.2 | 
2 
1 2 3 ~ 5 6 a 8 9 10 11 12 #13 #14 
Sample number 
Fig. 3. The effect of the freezing—thawing on density. 
= Before freezing After freezing 
3 
2.5 
-~ 2 
S 
& 
& 
z 15 


10 11 12 13 «14 


Fig. 4. The effect of the freezing—thawing on water absorption percentage. 


“al 
3 4 5 6 7 8 


9 
Sample number 
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= Before freezing After freezing 


2 3 ~ 5 6 7 8 9 10 11 12 13 14 


Sample number 
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Fig. 5. The effect of the freezing—-thawing on uniaxial compressive strength. 


= Before freezing After freezing 


2 3 A 5 6 7 8 9 10 11 12 13 #14 


Sample number 


100 
99.5 
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95 


DI (%) 


Fig. 6. The effect of the freezing—thawing on durability index. 


As can be seen in Table 3, the effects of freezing-thawing on the physical and mechanical 
characteristics are clear. After freezing-thawing, density, durability index, and uniaxial 
compressive strength are decreased and water absorption is increased. The percentage changes of 
studied parameters affected by freezing-thawing are given in Table 4 and Figure 7. 


As seen in Table 4, A5(Naein red granite) with a 13% reduction in the density highest, and A2 
(Azarshahr walnut travertine) and A3 (Piranshahr olive granite) have the lowest density 
reduction. Also, the average density of stones decreased by 7.02% due to freezing. It can be seen 
that Al (Khalkhal cream travertine) with 33%, the highest, and A4 (Mahabad gray granite) by 
3%, have the lowest increase water absorption, on average, water absorption percentage in the 
studied stones by freezing is 9.01%. 
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Table 4 
Physical and mechanical properties percentage changes of studied stone after freezing-thawing. 
Change (%) 
Sample ID 
Wa UCcS DI 
Ay -3 33 —13.86 —1.6 
Ay —2 8 —15.46 —1.21 
Ag =2 5 —6.63 —0.7 
Ay —8 3 —3.43 —1.5 
As —13 8 —4.78 —0.86 
Ag —6 9 —7.85 —0.94 
A = 16 fe =a 
Ag —10 6 —3.03 —1.16 
Ag —] 4 —8.57 —0.95 
Ajo —7.6 6 =3:15 —0.6 
Ay —7.1 7 —5.03 —1.22 
Ay 9 8 —13.28 —0.86 
Ay3 —11.4 4.5 —10.27 —0.91 
Ais —4.2 8.7 —8.72 —1.49 
--B--D(gr/em3) =- @—Wa(%) = 8 —UCS(MPa) —#®—DI(%) 
i 
30 \ 
25 \ 
20 s 
15 \ “A 
S 10 \ _ \ 
S ma -- ‘\ —_ o- 8. a 
ep 5 pe ee e_ = - ‘a? 
- 0 —_ a: St — | ee a ee ee si re 
. | f° 2 CR: Gy a eS oe Oe ks | 
i a ~% =m 4 a: . >> ents Tae ° Pus 
-10 7 ss. o” ‘-=- Laat = * 
7 ‘a’ ,~* & 
15 | "== 
-20 


Sample number 


Fig. 7. The percentage changes of studied parameters after freezing-thawing. 


Also, according to Table 4, A2 (Azarshahr walnut travertine) with 15.46%, the highest reduction 
in compressive strength, and A8 (Azarshahr yellow travertine) with 3.03%, have the lowest 
decrease in compressive strength in the effect of freezing. The mean decrease in the compressive 
strength of stones in the effect of freezing is 8.01%. Also, Al (Khalkhal cream travertine) with a 
1.57 drop in durability index has the highest, and A10 (Hamadan black granite) with a 0.6%, has 
the lowest decrease in durability index, in the effect of freezing the durability index drops into 
1.07% on average. 
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3. Determining the desirability degree of studied stones using PROMETHEE 


The PROMETHEE is a method for making multi-criteria decisions as an efficient solves method, 
using two words preferred and indifference to choose the best option. The method is applied in 
many fields of engineering such as operations research and dynamic management. This method 
has been well-attended because of its mathematical properties and its ease of use. According to 
the PROMETHEE method (Eq 1): 


Max (Min){ f(a), f,(@),--f, (a) | aeA} (1) 
A represents a set of decision options. 


The f,(a) shows (j) indicator value in option a, and j =1,2,...,k is the set of indicators or criteria 


that evaluates the options. Rating options by comparing the test is done on each indicator that 
this ranking is done in three steps as follows [42]: 


First step: Compare options using a pre-defined preference function with [0,+1] domain. For 
example, the preference function P, to examine options a and b and j indicator is as follows: 


P,(a,b) = P,| d, (a,b) | (2) 


So that: d, (a,b) = f(a) — f(b) represents the difference in the J indicator. Of course, for each fj 


indicator a weight (w}) also is considered. 


Second step: Calculating the overall preference of z(a,b) or each option. The higher z(a,b), the 
more preferable. The value of z(a,b)is calculated as follows: 


1(a,b) =) w;.p;(a.b) (3) 


k 
Where: > Ww, = 1, and J is weight indicator (w, ). 
j=l 
Third step: Calculating the overall preference of option a on other options. In this section output 
flow and also preferred the other choices on the option a is calculated as follows: 
for positive or output flow ratings: 


F 1 

O* (a) = —— > aa; x) (4) 
n—-1 xeA 

It shows how much priority an option a has over others. It shows the real power of option a. The 

@®*(a) largest amount means the best option. 

for negative input flow ratings: 


®O (a)= Ysa) (5) 


xeA 
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This shows that other options to what extent have priority over option a. In reality, this shows 
option a’s weaknesses. The smallest © (a) represents the best option. 

Minor ratings of options can be done by positive or negative flow PROMETHEE I. 

For complete rankings of all options the net flow ranking must be defined for each option 
PROMETHEE II: 


(a) = ©* (a) + © (a) (6) 


The resulting flow balance of positive and negative and net flow ranking higher, indicating a 
better option. For ranking and determining the desirability of studied stones by PROMETHEE 
the steps are to follow: 


Step 1. Forming decision matrix: This matrix indicates studied stones and their percentage 
change in density, water absorption, durable creep, and uniaxial compressive strength during 
different freezing—thawing cycles (Table 4). 


Step 2. Forming normalized decision matrix: For normalizing the data relating the value of 
each element is divided into the square root of the sum of squared elements in each column. The 
normalized decision is shown in Table 5. 


Table 5 
Normalized decision matrix. 
Sample ID D Wa UCS DI 

Ay 0.108 0.760 0.417 0.384 
A> 0.072 0.184 0.465 0.290 
A3 0.072 0.115 0.199 0.168 
Ag 0.287 0.069 0.103 0.359 
As 0.467 0.184 0.144 0.206 
Ao 0.215 0.207 0.236 0.226 
Ay 0.072 0.369 0.220 0.264 
Ag 0.359 0.138 0.091 0.276 
Ao 0.251 0.092 0.258 0.226 
Ajo 0.273 0.138 0.113 0.143 
Au 0.255 0.161 0.151 0.291 
Ap 0.323 0.184 0.400 0.205 
Ap 0.409 0.104 0.309 0.217 
Ais 0.151 0.200 0.262 0.356 


Step 3. Forming weighted normalized decision matrix: Each element in relating column in 
normalized matrix multiplied in the weight of every indicator. In this part the weight is assumed 
variable. Table 6 indicates the equal weight index weighted normalized decision matrix. 
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Table 6 
Weighted normalized decision matrix. 
Sample ID D Wa UCS DI 

Ay 0.027 0.190 0.104 0.096 
Ay 0.018 0.046 0.116 0.072 
Ag 0.018 0.029 0.050 0.042 
Ay 0.072 0.017 0.026 0.090 
A; 0.117 0.046 0.036 0.051 
Ag 0.054 0.052 0.059 0.056 
A; 0.018 0.092 0.055 0.066 
Ag 0.090 0.035 0.023 0.069 
Ag 0.063 0.023 0.064 0.057 
Axy 0.068 0.035 0.028 0.036 
Ay 0.064 0.040 0.038 0.073 
AY 0.081 0.046 0.100 0.051 
Ay 0.102 0.026 0.077 0.054 
Avs 0.038 0.050 0.066 0.089 


99 


Step 4. Calculating P,(a,b) and z(a,b): In this part after calculating differences of options, 


the amount P,(a,b) for each of the options about each indicator according to the first function 


would be achieved. In continue values for the options were calculated according to overall 


priority. The results of this review for option 1, according to Eqs 2 and 3, are shown in Table 7. 


Table 7 


Option | comparison with other options using the first z7(a, b) function with the value. 


D Wa UCS DI (a,b) 

(A;,A2) 0.009 0.144 — 0.012 0.023 0.75 
(Ay,A3) 0.009 ~=0.161 0.054 0.054 1 

(Ay,A4) —0.045 0.173 0.078 0.006 0.75 
(Ay,As) —0.090 0.144 0.068 0.044 0.75 
(A1,Ac6) —0.027 0.138 0.045 0.040 0.75 
(A1,A7) 0.009 0.098 0.049 0.030 1 

(Ay,As) —0.063 0.155 0.081 0.027 0.75 
(A1,Ao) —0.036 0.167 0.040 0.039 0.75 
(Ay,Ajo) —0.041 0.155 0.076 0.060 0.75 
(Ay,Aq1) —0.037 0.150 0.066 0.023 0.75 
(Ay, Aq) —0.054 0.144 0.004 0.045 0.75 
(Ay,A43) —0.075 0.164 0.027 0.042 0.75 
(Ay,Aq4) —0.011 0.140 0.039 0.007 0.75 


Step 5. Calculating minor ranking input and output flow: For selecting the best option, n—1 
options of A collection must be rejected. For this purpose, two priority flows are defined. 
Prioritizing output flow shows how much each choice has priority over other choices. 
Prioritizing input flow also shows how much of a particular option is superior to other options. 
So, the largest output and lowest input belong to the best option. Results of calculating the 
prioritized using Eqs 4 and 5 are shown in Table 8. 


Step 6. The net flow calculation and rating options: PROMETHEE I method calculates net 
flow options rating and represents each option’s superiority over other options. In this case, 
larger net flow represents the option’s superiority. The result of the calculation of net flow 
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prioritization using Eq 6 and ranking of studied dimension stones are shown in Table 8, Figures 8 
and 9. 


Table 8 
Input and output flow and net flow and final options ranking. 
Sample ID Output flow (positive) Input flow (negative) Net flow _ Ranking 


Ay 0.788 0.212 0.577 11 
Ad 0.558 0.365 0.192 9 
A3 0.173 0.788 — 0.615 1 
Ag 0.423 0.577 — 0.154 4 
As 0.500 0.462 0.038 7 
Ag 0.538 0.462 0.077 8 
Ay 0.481 0.481 0 6 
Ag 0.442 0.538 — 0.096 5 
Ag 0.404 0.596 — 0.192 3 
Ai 0.269 0.712 — 0.442 2 
Aun 0.519 0.481 0.038 7 
An 0.577 0.385 0.192 9 
Ay 0.538 0.462 0.077 8 
Aus 0.654 0.346 0.308 10 


4 — @ Output flow 
— @ — Input flow 


11 45 —*— Net flow 


Fig. 8. Input, output, and net flow of studied dimension stones. 


1 


8 
Fig. 9. Ranking of studied dimension stones. 
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As mentioned, indicators had the same importance. To review indicators’ importance in ranking 
in studied samples, sensitivity in indicators weight is considered. The results of the review 
showed in 13 states shown in Table 9. 


Table 9 
Ranking the options according to different weights assigned to the indicator. 


Indicators weight 


D 0.25 0.500 0.750 1 0.166 0.084 0 0.167 0.083 0 0.167 0.083 

Wa 0.25 0.167 0.083 0 0.05 0.750 1 0.166 0.084 0 0.187 0.083 
UCS 0.25 0.167 0.083 0 0.167 0.083 0 0.500 0.750 1 0.166 0.084 

DI 0.25 0.166 0.084 0 0.167 0.083 0 0.167 0.083 0 0.500 0.750 

ae Net flow 

Ai 0.577 0.205 —0.167 —0.538 0.719 0.858 1 0.666 0.757 0.846 0.718 0.859 

Ad 0.192 0.154 0.500 0.846 0.206 0.217 0.231 0.461 0.731 1 0.266 0.321 0.385 
A3 0.615 0.699 0.770 0.846 0.590 0.564 0.538 0.487 0.359 0.231 0.693 0.769 

Ay 0.154 0.024 0.206 0.385 0.436 0.718 -l 0.384 0.616 —0.846 0.180 0.512 0.846 
As 0.038 0.359 0.679 1 0.102 0.168 0.231 0.154 0.346 0.538 0.153 0.347 

Ae 0.077 0.025 0.129 0.231 0.282 0.487 0.692 0.076 0.078 0.077 0.026 0.128 

Ay 0 0.282 0.564 0.846 0.283 0.564 0.846 0.026 0.050 0.077 0.026 0.128 0.077 
Ag —0.096 0.166 0.430 0.692 0.167 0.236 0.308 0.397 0.699 -l 0.014 0.121 0.231 
Ag 0.192 0.154 0.115 0.077 0.410 0.628 0.846 0.051 0.089 0.231 0.154 0.115 

Ato —0.442 0.218 0.006 0.231 0.398 0.352 0.308 0.526 —0.609 0.692 0.628 0.815 

Au 0.038 0.051 0.065 0.077 0 0.038 0.077 0.102 0.244 0.385 0.205 0.372 0.538 
Av 0.192 0.308 0.422 0.538 0.205 0.218 0.231 0.359 0.526 0.692 0.103 0.397 

A 0.077 0.334 0.590 0.846 0.180 0.435 0.692 0.231 0.384 0.538 0.077 0.230 0.385 
Aus 0.308 0.077 —0.154 —0.385 0.385 0.461 0.538 0.333 0.359 0.385 0.436 0.564 0.692 


4. Particle Swarm Optimization 


In recent decades, artificial intelligence (AI) techniques, have played a significant role in solving 
complex problems in industry and science in theory and practice [46-53]. Numerous studies 
were conducted in a wide range of sciences that demonstrated the practical applications of these 
methods [54-59]. The Particle Swarm Optimization (PSO) is one of the most widely used and 
effective algorithms in solving various problems in conditions of uncertainty, which can be 
considered in the Swarm Intelligence (SI) which was first proposed by Kennedy and Eberhart 
(1995) [60]. This algorithm is based on the pattern of mass movement of fish in water or birds 
during migration to find its way to evolution and achieve the most optimal answer or answers. In 
this algorithm, starting with the optimization process, a set of particles (answers) are generated 
randomly. These particles move in set to get the best possible position in a D-dimensional space 
[61-63]. Each particle is represented by two vectors of position and velocity that represent the 
direction of its flight and velocity based on Figure 11. 


102 R. Mikaeil et al./ Journal of Soft Computing in Civil Engineering 6-1 (2022) 88-109 


Particle 
Oo 


New Position of Particle 


Particle Velocity 
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Fig. 11. New velocity and position vectors of a set of particles. 


Each particle, after moving at any moment, the velocity and position vectors are updated 
according to Eqs 7 and 8, respectively. Two components of the velocity vector are the personal 
best position (Pbest) and global best position (Gbest). 


ye = wv," +C 1,.( pbest, -X})+e, r,.(gbest — X*) 7) 
xe - x +Vi (8) 


Where V;*” and X\“*” represent the new velocity vector, and new position. Also, V,‘, and X} 
express the initial velocity vector and position of particles. The 7, and 7, are random numbers 


with uniform distribution with a range between 0 and 1. The w expresses the inertia weight to 
ensure convergence in the particle set. In addition, c,and c, are positive constants that called the 


individual learning factor and the social learning factor, respectively. Generally, the values of c, 
and c, are considered equal to 2, while Eq 9 must be met [64,65]. 


c+c¢, <4 (9) 


Position and velocity vectors of particles constantly being updated, and the optimization process 
is stopped when the best possible answer is achieved. It is worth mentioning that the PSO 
algorithm is used to train other algorithms such as artificial neural networks, neuro-fuzzy 
inference systems, and support machine vectors. Hence, in this study, the particle algorithm is 
used in training the Lloyd algorithm for clustering and data analysis, which is described in the 
following section. 


5. Clustering model and discussion 


As mentioned earlier, the PSO has a wide range of applications. Therefore, in this study, the PSO 
is employed to train Lloyd’s algorithm (k-means clustering). Clustering is one of the most 
effective methods in data analysis that the particle algorithm with the training of the Lloyd 
algorithm (one of the most popular clustering algorithms) increases the ability to classify analyze 
data to a desirable level [66]. Lloyd's algorithm is considered as a fitness function that is shown 
in Eq 10. 
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n . 
Obj.Function a min, d(x;,m p (10) 


Where d(xi,mj) represents the Euclidean distance of each member from each class’s center. j 
expresses the amount of clusters j=(1,2,3,...,k)[67—69]. 


For maximum algorithm efficiency in data analysis, the algorithm's control parameters must be 
determined that play an effective role in algorithm convergence. Since there are no specific 
relationships to determine these control parameters, a set of these parameters were selected based 
on previous studies [70-72]. Then, according to the data set and after trial and error, the most 
appropriate ones were selected including the number of clusters c=2, minimum acceptance 
precision of eL=0.00001, the initial population of 100, and maximum iteration of 150. The rock 
samples are classified into two classes based on the amount of physical and mechanical 
properties percentage changes after freezing-thawing. The results of this classification are 
demonstrated in Table 10 and compared with the results of ranking the rocks based on the 
permutation method. 


Table 10 
Comparison of the results of clustering by PSO in two classes with the results of ranking by 
PROMETHEE. 


PROMETHEE technique 
Sample ID Clustering by PSO 
Ranking Net flow 
Ai 2 11 0.577 
Ad 1 9 0.192 
A3 1 1 — 0.615 
Ag 1 4 — 0.154 
As 1 7 0.038 
Ag 1 8 0.077 
Ay 1 6 0 
Ag 1 5 — 0.096 
Ag 1 3 — 0.192 
Ato 1 2 — 0.442 
An 1 7 0.038 
Ay 1 9 0.192 
A 1 8 0.077 
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According to Table 3, the results of clustering by PSO algorithm showed that the rock sample Al 
was classified alone in one class and the other samples in another class together. Also, in a 
comparison with the results obtained from the PROMETHEE method, it is determined that the 
sample of Al, which is classified separately in the second class, has the lowest rank and has the 
highest amount of net flow compared to other samples. 


Finally, the following remarks can be concluded: 


The results of the experiments showed that after freezing density, uniaxial compressive strength, 
and durability index of stone are reduced and water absorption is increased. The reason for these 
changes can be analyzed in continuous periods of freezing and thawing that process of freezing 
and thawing sequence created and a series of large and small joints and cracks in the stone. 
However, if the joints and cracks appear on the surface small pieces of stone sometimes appear 
in the form of powder. If the cracks are within the volume of stone and the stone porosity is 
increased. The combination of these states reduces mass and increases the volume of sample 
stones and finally, characteristics such as density, strength, and durability of stone are decreased. 
Also, due to the formation of cracks in the stone, stone absorbs more water in free spaces and 
water absorption is increased. Studies on the rock-mechanical test of stone samples showed that 
AS5 (Naein red granite) with 13% decrease in density is highest and three walnut travertine, A3 
(Piranshahr Granat) and A7 (Tekab White travertine) with 2% have lowest density decrease 
during freezing. Also, in average stones density in freezing conditions decreased by 7.02%. 
Differences between the highest and lowest density percentage among studied stones were 11%. 
Al (Khalkhal cream travertine) with 33% increase and A4 (Mahabad gray granite) with 3%, in 
order, had the highest and lowest water absorption percentage increase. On average, the water 
absorption percentage showed a 9.01% increase in freezing conditions. Differences between the 
highest and lowest water absorption percentages among studied stones were 3%. A2 (Azarshahr 
walnut travertine) with 15.46% and A8 (Azarshahr yellow travertine) with 3.03 in order had 
highest and lowest compressive strength in freezing condition. The average reduction in stones 
compressive strength was 8.01. Also, Al (Khalkhal cream travertine) with 1.57% durability 
Index drop was highest and A10 (Hamadan black granite) with 0.6% had the lowest drop and in 
average durability index in freezing condition drops by 1.07. Then, using PROMETHEE multi- 
criteria decision-making with varied weights, the 14 samples were used to determine the degree 
of desirability of the stones. In most circumstances, the options A3, Al0, A8, and A4 have the 
maximum degree of freezing—thawing desire. Finally, classifying of laboratory data using PSO 
algorithm and comparison with the results of PROMETHEE ranking showed that the used 
methods can be considered as a suitable tool for assessing dimension stone according to 
resistance to freeze-thaw cycling to use in cold regions. 


6. Conclusions 


Freezing is one of the most effective environmental factors to the physical and mechanical 
properties of dimension stones. Since freezing is a destructive agent, thus causes the undesirable 
stone conditions and reduces quality and its efficiency. In this study, freezing effects on the 
physical and mechanical dimension stones were evaluated and the desirability degree of these 
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stones was ranked according to four percentage changes of important physical and mechanical 
such as density, water absorption, uniaxial compressive strength, durability index. For this 
purposel4 rock types of the most widely used dimension stones in cold regions were collected 
and transferred to the laboratory. To analyze the actual conditions of freezing in the cold areas 
that face freezing in winter, stones were studied at —20°C. Determining the degree of desirability 
studied stones using PROMETHEE multi-criteria decision-making with different weights for the 
14 cases studied and evaluated. The results of this study showed that in most cases the option A3 
(Piranshahr Granat), A10 (Hamadan black granite), A8 (Azarshahr yellow travertine), and A4 
(Mahabad gray granite) have the highest degree of desirability in terms of freezing—thawing. 
Thus, sample stones of A3, Al0, A8, and A4 are the most desirable stones in cold climates 
among studied stones. Then, the results of the laboratory were classified by PSO algorithm and 
compared with the results of PROMETHEE. Consequently, it could be recognized as the lowest 
rank in a separate class. This study indicated that the proposed approach can be applied as a 
reliable tool for evaluating freezing phenomena in rock mechanics engineering. 
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